arXiv:hep-ph/0112115vl 7 Dec 2001 


Radiative dileptonic decays of B-meson in the general two 

Higgs doublet model 

G. Erkol and G. Turan 

Physics Department, Middle East Technical University, 06531, Ankara, Turkey 
E-mail; gurerk@newton.physics.metu.edu.tr 
E-mail: gsevgur@metu.edu.tr 


We investigate the exclusive B decay in the general two Higgs 

doublet model (model III) including the neutral Higgs boson effects with an em¬ 
phasis on possible CP-violating effects. Eor this decay, we analyse the dependen¬ 
cies of the forward-backward asymmetry of the lepton pair, Aps, CP-violating 
asymmetry, Acp, and the CP-violating asymmetry in forward-backward asym¬ 
metry, Acp{Apb), on the model parameters and also on the neutral Higgs boson 
effects. We have found that App 10“^, 10“^, Aqp ^ 10“^, 10“^ and 
Acp{Apb) ^ 10“^ , 10“^ depending on the relative magnitude of the Yukawa 
couplings ^ and in the model III. We also observe that these physical 
quantities are sensitive to the model parameters and neutral Higgs boson effects 
are quite sizable for some values of the coupling 

PACS number(s): 12.60.Pr, 13.20.He 


1. Introduction 

It has been realized for a long time that rare B-meson decays induced by the 
flavor-changing neutral currents (FCNC) are the most promising field for obtain¬ 
ing information about the fundamental parameters of the standard model (SM), like 
the elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, the leptonic de¬ 
cay constants etc., and testing the SM predictions at loop level. At the same time 
rare decays can also serve as a good probe for establishing new physics beyond 
the SM, such as the two Higgs doublet model (2HDM), minimal supersymmet¬ 
ric extension of the SM (MSSM) [1], etc., since the contributions from these new 
models and the SM arise at the same order in perturbation theory. The observation 
of radiative penguin mediated processes, in both the exclusive B K*'y [2] and 
inclusive B Xgj [3] channels, stimulated the investigation of the radiative rare 
B meson decays with a new momentum. Among these rare decays, B —> ji+r 
{i = e, fi, t) have received a special interest due to their relative cleanliness and 
sensitivity to new physics. They have been investigated in the framework of the 
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SM for light and heavy lepton modes in refs. [4]-[6]. The new physies effeets in 
these deeays have also been studied in some models, like MSSM [7] and the two 
Higgs doublet model [8, 9, 10]. 

In this work, we study the radiative B ^ j deeay in the general two- 
Higgs doublet model (model III) ineluding the neutral Higgs effeets. The 2HDM 
is the minimal extension of the SM, which consists of adding a second doublet 
to the Higgs sector. In this model, the Yukawa Lagrangian responsible for the 
interaction of quarks and leptons with gauge bosons opens up the possibility of 
having tree-level FCNC, which are forbidden in the SM and model I and II types 
of the 2HDM. This brings new parameters, i.e., Yukawa couplings, into the theory. 

H ^ 7 decay is induced by the pure-leptonic decay B , which is 

free from the helicity suppression, in contrast to the channels with light leptons, but 
quite hard to detect experimentally due to low efficiency. In H ^ 'y decay, 
helicity suppression is overcome by the photon emission in addition to the lepton 
pair. For this reason, it is expected for H ^ 7 £~^£~ decay to have a large branch¬ 
ing ratio and this makes its investigation interesting. Another reason that motivates 
to study H ^ 7 £~^£~ process is that it receives additional contributions from the 
neutral Higgs boson (NHB) exchanges in the 2HDM. Since NHB contributions 
are proportional to either the lepton mass or the corresponding Yukawa coupling, 
they are negligible for B -f£+£- decays with light leptons, but we could expect 
significant contributions for £ = t. Indeed, the investigation of B ^ 
decay in model I and II types of the 2HDM in [8], and in MSSM in [7], including 
NHB effects report that the contribution from exchanging neutral Higgs bosons 
may be quite sizable for large values of tan (5, which is already favored by LEP 
experiments [11]. 

We investigated the B —> -f£+£- decay for £ = r in the model III type of 
the 2HDM in a previous paper [10]. In this work we extend this study with an 
emphasis on possible CP violation effects. The CP asymmetry is of great interest 
in high energy physics especially since its origin is still unclear. In the SM, the 
source of CP violation is the complex CKM matrix elements, which can explain 
all the existing data on CP violation. However, for example, to explain the matter- 
antimatter asymmetry observed in the universe today one needs additional sources 
of CP violating effects, which has motivated to search new models beyond the SM. 
In model III type of the 2HDM, the complex Yukawa couplings provide a possible 
source of CP violation. Indeed, it was reported [12]-[14] that a measurable CP 
aysmmetry was obtained due to this new phase in the model III. 

The paper is organized as follows: In section 2, we first present the lead¬ 
ing order (LO) QCD corrected effective Hamiltonian for the quark level process 
b 7 including the NHB exchanges. Then we give the corresponding 

matrix element for the exclusive B ^£'^£~ decay. Next, we calculate the 
the forward-backward asymmetry of the lepton pair, ApB, CP-violating asym¬ 
metry, Acp, and the CP-violating asymmetry in forward-backward asymmetry. 
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Acp{Afb), as functions of the model parameters. Section 3 is devoted to the 
numerical analysis of these physical quantities with respect to the CP parameter 
sin 6, Yukawa couplings and and the mass ratio m/jO/m^o and to the 
discussion of our results. In Appendix A, we give a brief summary about the gen¬ 
eral 2HDM (model III). The operators and the corresponding Wilson coeffients 
appearing in the effective Hamiltonian are given in Appendices B and C, respec¬ 
tively. Finally, some parametrizations used in the text may be found in Appendix 
D. 


2. The exclusive S —> 7 £+£ decay 

The exclusive 5^7 decay is induced by the inclusive b ^ sj 
one. Therefore, we start with the QCD corrected effective Hamiltonian for the 
related quark level process b ^ s , which is obtained by integrating out heavy 
particles in the SM and in the 2HDM [15, 16] 


n 


-4Gf 


VtbVt 


ts 


10 10 'I 


i=l 


i=l 


( 1 ) 


where Oi are current-current {i = 1 , 2 ), penguin {i = 1 ,.., 6 ), magnetic penguin 
(i = 7, 8 ) and semileptonic (i = 9,10) operators . The additional operators 
Qi,{i = 1,.., 10) are due to the NHB exchange diagrams, which give considerable 
contributions in the case that the lepton pair is [15]. Ci{ii) and CQ^{|J,) are 
Wilson coefficients renormalized at the scale /U. All these operators and the Wilson 
coefficients, together with their initial values calculated at = mw in the SM and 
also the additional coefficients coming from the new Higgs scalars are presented 
in Appendices B and C. 

Neglecting the strange quark mass, the effective Hamiltonian (1) leads to the 
following matrix element for 6 ^ s £^£~, 


M 


^VtbVZ}^Cl^f{s^^PLb) + CMs^^Pib) h,j,£ 
2C7^{siaf,^q„PRb)£-f^£ + CQ^{sPRb)££ + CQ^{sPRb)£-f5£\ , (2) 


where Pl,r = (l=F 75 )/ 2 ,qis the momentum transfer and Vij’s are the corre¬ 
sponding elements of the CKM matrix. 

In order to obtain the matrix element for 6 ^ 57 £^£~ decay, a photon line 
should be attached to any charged internal or external line. As pointed out before 
[5, 6 ], contributions coming from the release of the free photon from any charged 
internal line will be suppressed by a factor of m^/M^ and we neglect them in 
the following analysis. When a photon is released from the initial quark lines 
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it contributes to the so-called ’’structure dependent” (SD) part of the amplitude. 
Using the expressions 


(7(^)|s7m(1=F 75)^1 


m. 


{e^.XaS*’^q^k^g{q^)±i[s*^^{kq) - is*q)k^^]fiq^)} 


= —{ef^apae^kf^qa giip^) =F i ^l.{kq) - {e*k)q^ , 


rriB ^ 

{l\s{l + l5)b\B) = 0, 

the SD part of the amplitude can be written as 

cxGp 


( 3 ) 


Msn = 


2x/2 


-ytbVts-^ S Mei,uape*''q°‘k^ + iMi£l{kq) - is*q)kA 

TT I ^ ^ 


+ 


irid [Bie^,^p£*'^q^k^ + iB 2 (el{kq) - {e*q)k^) 


( 4 ) 


where e* and are the four vector polarization and four momentum of the photon, 
respectively, and pb is the momentum of the B meson. In Eq. (4), Ai, A 2 , Bi 
and B 2 are functions of the Wilson coefficients and the form factors, and they are 
given in Appendix D. 

We note that the neutral Higgs exchange interactions do not contribute to the 
structure dependent part of the matrix element M.sd- However, the situation is dif¬ 
ferent for the so-called ’’internal Bremsstrahlung” (IB) contribution, M.ib, which 
arises when a photon is radiated from one of the final f- lepfons. Using fhe expres¬ 
sions 


(0|s7/^75^|^) 
{0\saf,u{l ± j5)b\B) 

(0|s756|H) 

m\B) 


- ifsPBn , 

0 , 


ruh 


0 , 


( 5 ) 


and fhe conservation of fhe vector currenf, IB pari of fhe malrix elemenf is found 
as [8] 


Mib = 


aGp 

2x/2 


TT 


VMefBi Fi 


f i>B _ i’B^' 
y 2piA: 2-p2k 


175^ 


2pik 2p2k 


+ 2m£ 


1 1 

-h 


2pik 2p2k 




+ Fii 


( 6 ) 
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where F and Fi are funetions of form faetors and the Wilson eoeffieients Cq^ and 
Cq 2 due to the NHB effeets and their explieit forms ean be found in Appendix D. 
Finally, the total matrix element for the B deeay is obtained as a sum 

of the M.SD and Mib terms, Ad = M.sd + -Mib- 

Now, we will ealeulate the forward-baekward asymmetry, Aps, for the lepton 
pair, CP-violating asymmetry, Acp, and CP violating asymmetry in the forward- 
baekward asymmetry, Acp{Afb) for the proeess under eonsideration. All these 
measurable physieal quantities ean provide a great deal of elues to test the theoret- 
ieal models used. We first give the definitions of Afb{x) and Acp- 


Afb{x) 


Acp 


T{B je+i-) - T{B ji+i-) 
r{B ji+i-) + r(.B ^ -fi+i-) 


(V) 

( 8 ) 


where 2 : = cos 6, 0 is the angle between the momentum of the B-meson and that 
of i~ and X = 2Ey/mB is the dimensionless photon energy. In Eq.(7), 


d'^r 


IS 


the double differential decay rate and in the center of mass (CM) frame of the 
dileptons it is given by 


d^T 

dxdz 


aGp 


2^/2 


VtbVt 


TT 


ts 


a TT 


-n — niB X V 

{2Ff 4 


m 


B ^2 


32 


((1 + Z^)(l-x- 4r))(|Ai|2 + 1^2^ + \Bi\^ + W) 


+ 8r(|Ai|2 + lAsI^) + 42^(1 - x)(l - x - 4r)Re(A2Si* + AiB^) 
(x-l) 


+ fpryii- 


vxzRe{B2F* — BiF^ 


((z^ — l)(x — 1) + 

+ (1 - 4r - 2^(1 - X - 4r))Re(A2Fi) - xRe(AiF*) 
(1-x) 


+ fB- 


\F\ ((—2 + 4x — 3x"^ + x'^ 


' x^{{z^ — l)(x — 1) + Arz^)^ 

+ 8r(l — x)){z^ “ 1) + 4rx^z^^ + \Fi\^ ^^32r^(x — 1) + 4r(4 — 6x + 2x^) 


— 2 + 4x — 3x^ + x^^ {z^ “ 1) + x^z^'j I . 


where v = 


\J^~ 1 ^ ''^hh r = 


Integrating over the angle variable, we 


( 9 ) 
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find the forward backward asymmetry ApB follows, 


ApB = —J \j{x — l)(x — 1 + 4r)Re(Aiyl2 “ B1B2) 

I nf] — 1 \ At 

- " ^-l + 4r U Y^Re((4l, - B 2 )F- 


- (^1 -Bi)F;)| j J dxD{x), 


where 


D{x) = (|Ai|2 + |A 2 |")(l + 2r-x) + (|Bi|2 + | 52 p)(l- 4 r-x) 

— fsiTitx 2v{\ — x)Re(^ 2 -^r) + ln | ^ ^(x — 4 r)Re(A 2 -^r) — xRe(^iF*)^ 

- 2/1 [„il^ ( xr . a - 4 .) IF,, j F .n^ ((: F I - i - I) ,F,^ 

+ ('(l_4r)-MLl^±^_|')|F,|^)], (11, 


We note that in these integrals the Dalitz boundary for the dimensionless photon 
energy x is taken as 

5 <x<l- ^ , (12) 

m% 

since \M.ib\^ term has infrared singularity due to the emission of soft photon. 
In order to obtain a finite result, we follow the approach described in ref. [ 6 ] and 
impose a cut on the photon energy, i.e., we require > 50 MeV, which cor¬ 
responds to detecting only hard photons experimentally. This cut requires that 
> bruBj^. with b = 0 . 01 . 

For B 7 decay, Acp almost vanishes in the SM due to the unitar- 
ity of CKM matrix together with the smallness of VubV*g. However, in model III 
complex Yukawa couplings provide a new source of CP violation. In our calcu¬ 
lations, we choose '^hat , Cq^ and Cq^ are the Wilson 

coefficients that contain CP violating terms. Using Eq.( 8 ), we calculate Aqp as 


Acp 


J dxT{x) 

f dx {D{x) + Dcp{x)) ’ 


(13) 
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where 


T{x) = m% irn{^N,bb) \^x{l + 2r - x)A^^^ 


-Vb — 

mb 


(2v(l — x) + (x — 4r)ln— A^^FP^ — xln— A?'^\ , 
\ 1 — V J 1 — V JJ 


(14) 


and Dcp{x) is the CP eonjugate of D{x) whieh is defined as 

Dcp{x) = D{x) ^ . (15) 

The explieit form of the funetions ^(ij’ ^^l-Cl^) are given in Appendix 

D. 

Finally, we eonsider the CP violating asymmetry m App, Acp{Afb), whieh 
is an important measurable quantity that may provide information about the model 
used. It is defined as 


Acp{Afb) 


Afb —Afb 
Afb+Afb 


(16) 


Here, Afb is given by Eq.(lO) and Afb is obtained by the replaeement ^^bb 
i^N,bb)* “ ^FB ■ 


3. Numerical analysis and discussion 

We present here our numerieal results only for £ = r ehannel, but they ean 
easily be applied to the £ = fi ease. The input parameters we used in our numerieal 
analysis are as follows: 

mB = 5.28 GeV , mb = 4.8 GeV , me = 1.4 GeV , mr = 1.78 GeV , 
mjjo = 150 GeV , m^o = 70 GeV , m^o = 80 GeV , mjj± = 400 GeV , 

\VtbVZ\ = 0.045, = 129, Gp = 1.17 x 10“^ , tb = 1.64 x 10“^2 s , 

= 4.229 , Cio = -4.659. (17) 

Here we note that the value of the Wilson eoeffieient Gg^^ above eorresponds 
to only the short-distanee eontributions. Gg^^ also reeeives long-distanee (LD) 
eontributions assoeiated with the real cc intermediate states (See Appendix C for 
the details of LD eontributions). There are five possible resonanees in the cc system 
that ean eontribute to the H ^ 7 deeay and to ealeulate their eontributions 
we need to divide the integration region for x into two parts: (5 < x < 1 — ((m ^2 + 



0.02)/mB)^ and 1 — ((m ^2 “ 0.02)/ms)^ < x < 1 — {2mr/mB)‘^, where 
m ^2 = 3.686 GeV is the mass of the seeond resonanee. 

For the values of the form faetors g, /, gi and fi, we have used the results of 
ref. [17] and [18], and represent their dependeneies in terms of two parameters 
F{0) and mp ss 


n<f) = 



(18) 


where the values F{0) and mp for the B ^ ^ are listed in Table 1. 

There are many free parameters in the general 2HDM, sueh as masses of the 
eharged and neutral Higgs bosons and eomplex Yukawa eouplings, where 

i,j are quark flavor indiees. There are also some experimental results that one 
ean use to restriet these new parameters. In this eontext, the stronger restrietion 
eomes from the analysis of the AF = 2 deeays with F = K,Bd,D mesons, the 
p parameter and the B ^ Xg'y deeay. 

The eontributions to the Wilson eoeffieient Cj from the neutral Higgs bosons 
and are given by [19] 


Cj{mw) 


i=d,s,b 


Qi 

Sniimb ’ 


(19) 


where H = , A^, and m* and Qi are the masses and eharges of the down quarks, 

respeetively. Eq. (19) shows that the neutral Higgs bosons ean give a large eontri- 
bution to the Wilson eoeffieient C 7 and this eontradiets with the CLEO data [3] 

BR{B ^ Xs-f) = (3.15 ±0.35 ±0.32) X 10“^ (20) 

However, assuming that the eouplings ig with i = d,s,b and small 

enough to reaeh the eonditions (N,bb^N,is ^ CN,bd^N,ds ^ ^ (1^], together 
with the eonstraints from AF = 2 mixing [20] and the p parameter [21], we obtain 
the eonditions 

iN,tc ^ ) 

^N,ib ) ^N,ij ^ 0i,j=d,S, (21) 

SO that we only take into aeeount ^ and In our work, we ehoose ^ 

as real and complex, namely i^bb — the ^ sinee 

it eontrols the eontributions due to the NHB effeets, we leave ^ as a free pa¬ 
rameter to investigate the size of NHB effeet on the measurable quantities of the 
B ^T^T~ deeay. 
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In our numerical calculations, we further adopted the constraint on the Wilson 
coefficient , 0.257 < < 0.439 [20] due to the CLEO measurement 

(Eq.(20)) and the redefinition 



The above constraint on the restricts the allowed regions of the measur¬ 
able quantities of the exclusive B decay, namely, Afb, Acp and 

Acp{Afb) and these regions are represented by the ones between the solid curves 
for > 0 and the dashed curves for < 0 throughout the graphs in Eigs. 
1-18. 

In Eig. 1, we plot sin 6* dependence of Afb without NHB effects for the case 
of ratio < 1. We see that in model III without NHB effects, | 

is smaller compared to its value in the SM (0.183), represented by the dashed 
straight line, for > 0, but it can be enhanced up to 7% with increasing 

sin0. Eor < 0, is not much sensitive to sin0, but its value can be 
slightly greater than the SM prediction. Including the NHB effects to Afb { Fig- 
2) reduces its magnitude 30% of its value without NHB effects for > 0, 
while for < 0, yliTB is almost the same as the SM value. 

In Eig. 3 (4), we present Afb ^ function of (m/jo/m^o) for = 
10 m,-, sin0 = 0.5 and |rib| < 1. Afb is at the order of magnitude 10“^ and 
increases with the increasing values of both and Eor > 0, 

Afb stands less than the SM prediction and for < 0, model III prediction 
can reach the SM one for large values of and m^o / nij^o. 

We have also calculated the sin 6* ( and mifi/rriAo) dependence of Afb 
for the case of ratio rtb > 1- We have found that in this case Afb one (two) 
order(s) of magnitude smaller than its value for \rtb\ < 1 case, and including NHB 
effects reduces this value even one more order of magnitude. Therefore we do not 
present these graphs here. 

Eig. 5 represents sin0 dependence of Acp without NHB effects for the case 
of ratio < 1. It is at the order of magnitude 10“^ and increases with sin 6 . Eor 
< 0, Acp can have both signs, while for > 0 its sign does not change 
in the restricted region. Including the NHB effects (Eig. 6) reduces I^cpI without 
NHB effects almost by 60% for > 0. However, for < 0, it is possible 
to enhance it by up to 35%. 

Eigs. 8 and 7 are devoted to sin0 dependence of Aqp for rtb > 1 with and 
without NHB effects, respectively. Without NHB effects, |Acp| is at the order of 
magnitude 10“^ and including the NHB effects can enhance it up to two orders 
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of magnitude, i.e., it becomes |^cp| ~ 10 We further note that the restricted 
region for > 0 (solid lines ) and < 0 (dashed lines) are now larger but 
they almost coincide. 

Fig. 9 (10) shows Acp as a function of for = 40(0.l)m;, and 
sin6* = 0.5 for \rtb\ < 1 (rtb > 1)- We see that Acp is sensitive to the parameter 

T-T- and it decreases (increases) with the increasing values of for > 0 
< 0) when |rtf,| < 1. When rtb > 1? restricted regions for > 0 (solid 
lines ) and < 0 (dashed lines) almost coincide and | Acp| can take one orders 
of magnitude larger values compared to the case where |rtf,| < 1. 

In Fig. 11 (12), we plot the dependence of Acp on for = 

10(l)mT-, = 40(0.1)mfe and sin0 = 0.5 for \rtb\ < 1 {rtb > 1)- For 

\rtb\ < 1 and > 0, Acp is sensitive to the ratio mf^o/mjfi and increases 
with the increasing values of m/jO /m^o. However, for < 0, dependence of 
Acp on the ratio is weak but, the restricted region is larger this time. As seen 
from Fig. 12, when rtb > 1^ I^cpI can take one orders of magnitude larger values 
compared to the case where |rtf,| < 1. 

Finally, we present our results about the CP violating asymmetry in App, 
Acp{Afb) in a series of figures. Figs.13-18. Fig. 13 (14) shows Acp{Apb) 
as a function of sin0 for = 10(1) mr, ^Nbb — 40(0.l)mb for \rtb\ < 1 
{rtb > 1)- For I rtf, I < 1, Acp{Afb) is at the order of magnitude 10“^ and it does 
not change sign in the restricted region for > 0, while it can have both signs 
for < 0. For rtf, > 1 , Acp{Afb) can reach 5% for the intermediate values 
of sin 6 and restricted regions for > 0 (solid lines ) and < 0 (dashed 
lines) almost coincide. 

We can see from Figs.15 and 16 that Acp{Afb) is sensitive to the parameter 
CnTT’ especially for its small values. Acp{Afb) is a decreasing (increasing) 
function of for |rtf,| < 1 {rtb > 1) and reaches 1(6)% for = 1(50). 

Fig. 17 (18) is devoted to the ratio dependence of Acp{Afb) for 

= 10(1) m^, i^ tib = 40(0.1)mf, and sin^ = 0.5 for |rtf,| < 1 {rtb > !)■ 
It is seen from Fig. 17 that for |rtf,| < 1, Acp{Afb) is sensitive to the mass 
ratio m/jO / rrij^o and it is increasing with the increasing values of m/jO / rrij^o for 
> 0, while for < 0, its dependence on the mass ratio is weak. For 
rtb > 1> \Acp{Afb)\ is one order of magnitude larger than its value for |rff,| < 1. 

In conclusion, we have investigated the physical quantities App, Acp and 
Acp{Afb) for the exlusive B j decay in the general 2HDM including 
the NHB effects. From the results we have obtained we conclude that experimental 
investigation of these quantities may be very useful for testing the new physics 
effects beyond the SM and also the sign of . 
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F(0) 

ap 

9 

IGeV 

5.6 Gel/ 

f 

0.8 GeV 

6.5 Gel/ 

91 

3.74 Gel/2 

6.4 Gel/ 

fi 

0.68 Gel/2 

5.5 Gel/ 


Table 1. B meson decay form factors in the light-cone QCD sum rule. 


Appendix A 

Model Description 


The 2HDM is the minimal extension of the SM, which consists of adding a second 
doublet to the Higgs sector. In this model, there are one charged Higgs scalar, two 
neutral Higgs scalars and one neutral Higgs pseudoscalar. The general Yukawa 
Lagrangian, which is responsible for the interactions of quarks with gauge bosons, 
can be written as 

= TlijQiL(i>lUjR +r]fjQiL(l>lDjR +QiL(j)2UjR +iijQiL(l)2DjR +kh.1), 

where i, j are family indices of quarks , L and R denote chiral projections L{R) = 
1/2(1 =F 75), (j)m for m = 1, 2, are the two scalar doublets, QiR are quark doublets, 
UjR, DjR are the corresponding quark singlets, and are the matrices 
of the Yukawa couplings. The Yukawa Lagrangian in Eq. (A.l) opens up the 
possibility that there appear tree-level FCNC, which are forbidden in the SM and 
model I and model II types of the 2HDM. However, tree-level FCNC are permitted 
in the general 2HDM, and this type of 2HDM is referred to as model III in the 
literature. 

In this model, it is possible to choose (pi and (j )2 in the following form 





0 

v + H^ 


+ 





V2H+ 
Hi + iH2 


,(A.2) 


with the vacuum expectation values, 

<</>!>= ^ ( ) ; < (/>2 >= 0 . (A.3) 

With this choice, the SM particles can be collected in the first doublet and the new 
particles in the second one. Further, we take Hi, H 2 as the mass eigenstates hP, 
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AP respectively. Note that, at tree level, there is no mixing among CP even neutral 
Higgs bosons, namely the SM one, and beyond, h^. 

The part which produces FCNC at tree level is 

^Y,FC = ^fj^QiL(j)2UjR + ^ijQiL4'2DjR + ^klhL4>2ElR + h.C. . (A.4) 

In Eq.(A.4), the couplings for the flavor-changing charged interactions are 

^ch — ^neutral ^CKM t 

^ch ^CKM ^neutral ; (A-5) 

where ^^eutrai defined by the expression 

<a.6) 

and Cneutrai denoted as Here the charged couplings are the linear com¬ 

binations of neutral couplings multiplied by Vckm matrix elements (see [20] for 
details). 


Appendix B 

The operator basis 

The operator basis in the general 2HDM (model III) for our process is [15, 16, 22, 

23] 

Ol = {sLalk-^L0){cL0l^bLa) , O 2 = (sLa7/iCL«)) 

O 3 = isLalk^La) {qL0l^qL0) , O 4 = {sLalk^L0) Y i^L0l^qLa) , 

q=u,d,s,c,b q=u,d,s,c,b 

O 5 = (sLalt^bLa) Y , ^6 = {sLaliibL 0 ) Y , 

q=u,d,s,c,b q=u,d,s,c,b 

O7 = Y^SadiMuirribR + msL)baiF^'' , Os = Y^SaT^0crf,u{mbR + msL)bpg°-> 
O9 = -:^^{sLaliibLa){h^t) , Oio = r^^isLalb^bRaWl^ld) , 

Qi = -^{si bi) m , Q 2 = bi) {hd) , 

Q^ = -^^^LbR) Y (^l 4 ) > <34 = II i^Rqi) ^ 


q=u^d,s,c^b 


q=u,d,s,c,b 


<35 = H (^L<1r) ^ = ^^i^lbR), Y 

q=u,d^s,c,b q=u,d,s,c,b 
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2 2 
y /-a ^uv i,a \ sr^ ^ /o _ y 


Qi = b%) J2 (^L > <38 = Y^i^L b%} J2 (4 ^ m ^ 4 ) > 

q=u,d,s,c,b q=u,d,s,c,b 

Q9 = r^islay’^b^jY Y1 ^ Qw = r^^isla^’^b^j^) ^ {fR^yi^yql) 


q=u,d^s,c,b 


lOvr^ 


q=u,d,s,c,b 


where a and (3 are SU (3) eolour indiees and and Q^'’' are the field strength 
tensors of the eleetromagnetie and strong interaetions, respeetively. Note that there 
are also flipped ehirality partners of these operators, whieh ean be obtained by 
interehanging L and R in the basis given above in model III. However, we do not 
present them here sinee eorresponding Wilson eoeffieients are negligible. 


Appendix C 

The Initial values of the Wilson coefficients. 


The initial values of the Wilson eoeffieients for the relevant proeess in the SM are 

[ 22 ] 


Ci^imw) 

Cj^{mw) 

Ci^{mw) 

Ci^imw) 

Cf^{mw) 


0 , 

1 , 

3xj — 2 x1 
4(xt - 1)4 
3x1 


Inxt + 


—Sxf — 5x1 + 7xt 


4{xt - 1)4 

1 


Inxt + 


24(xi - 1)3 
— xf + 5x1 + 


8 {xt - 1)3 


B{xt) + 


1 — 4 sin^ 6 w 


siv?9w 

-^^{B{xt) - C{xt)) 


-C'(xt) - D{xt) + - , 

sm 6 w 9 


sin^ 6*vy 

0 i = l,..,10 


(C.l) 


and for the additional part due to eharged Higgs bosons are 

Ci,...Q{mw) = 0, 

C^{mw) = Y^Fiiyt) + XY F2{yt) , 
C^{mw) = Y^Gi{yt) + XY G2{yt) , 
Gl^{mw) = Y^HYvt) , 

C'io("^'Vl^) = Y"^ Li{yt) , 


(C.2) 
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where 


1 


^ — ( CN,bb + iN,sbT^ 

mb \ ’ Vtb 

1 f -TJ ^JJ V*g 


.D Vts 


mt 




V, 


ts 


(C.3) 


The NHB effeets bring new operators and the eorresponding Wilson eoeffieients 
read as [24] 




'Q 2 


r^A°f/tU ^2^ _ ^N,TTi^N,tt) (n, 

) - 327r2m^o 


^N,TTi^N,tt) ^byt{&5{yt)zA - Qi{zA,yt)) 

327r‘^m\omt<di{zA,yt)&5{yt) 

D f cU \2cD -I 


bb , 


{zA,yt)&iizA,yt)&5iyt 


{yt{&i{zA,yt) - &5{yt){xy + za)) - 2ei{zA,yt)&5{yt)H ^ I ]) 

&i(ZA,yt) 


^ao( 7 U ^ _ 9 ^^N,TAN,ttbnbXt ( 2 _ xyxt + 2 za _ , . ZA^^ixt) . 

Q 2 N,tt MTx‘^m\omt ^©5(3;*) Qi{zA,xt) 0 i(zA,a:t) 

{x - l)xt{yt/zA - 1) - (1 + x)yt) 


- xyxtyti 


(06 - {x - y){xt - yt)){Q^{zA) + ix- y){xt - yt)zA) 
x{yt + xt{l - yt/zA)) - 

0603(^a) ’) 


cD 

riA^rFD \ ^ ^N,TT^N,bb I 

- 64^2^2^ I 


1 - 


xht + 2 y(x - l)xtyt - ZA{xt + ©e) 


03(^;4) 


Xti'^-yt/zA) ^ zaQ& 


06 


Q2 {za,x] 


^ 2 \ _ y^N,tt) 'bbibmr ( xt{l - 2 y)yt ^ (-1 + 2 cos^ %)(-! +x + y)?/^ 

) - 64^2^^^^2 I + 


&5{yt 


cos^ 9w&4{yt) 


+ 


ZHi&i{zH,yt)xyt + cos^ 0w {-2x‘^{-l + xt)yyt + xxtyyt - Q?,zh))\ 


cos2 6wQi{zH,yt)&7 


} 


r<H°(cU \ _ 9 ^N,tt^N,bb''^'r ( (-1 + 2cos2 Ow) yt xtyt , xtyt{xy - zh) 

^Qi \ZN,tt) — 


+ 


+ 


647r‘^mjjomt \ cos^ 04(yt) &5{yt) &i{zH,yt) 

{-1+ 2 cos'^ 9w)ytZH „ 1 \ Q5(.yt)zH 


COS^ 0 H /07 


— 2 xt In 


QiizH,yt) 


(C.4) 
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( 4\ _ 9 'ni^fTtlT-Xt 

^ 128 tt‘^ rnjjoTnf 


- 1 + 


(-1 + 2x)xt 

Q5{xt) + y(i -xt) 


2xti-l + {2 + xt)y) 

&5{xt) 


4cos^6'h/(- 1+ x + y) + xt(x + 2/) xt{x{xt{y - 2 zh) - 4:Zh) +‘2zh) 


COS^ 6*14/04(xt) 


&l{zH,Xt) 


+ 


yt{{-l + x)xtZH + cos^ Owii^x - y)zH + xt{2y{x - 1) - zh{ 2 - 3x - y)))) 


cos2 9w{&3{zh) + x{xt - yt)zH) 


+2 (xt In 


&5{xt)zH 

&l{zH,Xt) 


+ In 


x{yt - xt)zH - QsjzH) 

i&bixt) + 2/(1 -xt)ytZH 




327r^m^o 


^N,TTi^N,tt) - byt - 1) + 05 (yt)( 2 x - l)zh) 




^N^rri^N^ttf ( {Qb{yt)zh{yt “ l)(x + y- 1 ) - 01 (^/i, yt) (05 (//t) +yt) 


32 vr 2^2 


- 21n 


\ 

Zh^bjyt) 

Ql{Zh) 


Qi{zh)Qb{yt 


r-<h° \ _ 9^^N,TT^N,tt^bXt f 2{—l + {2 + xt)y) xt{x — l){yt — Zh) ^ 

w^) ^ 


Zh&djxt) 

&l{Zh,Xt) 


+ 


x{xt{y - 2zh) - 4zfe) + 2zh 

^l{Zh,Xt) 


(1 + x)ytZh 


xyxtyt + Zh{{x - y){xt - yt) - 06) 


€>9 + ytZh{{x - y){xt - yt) - 06)(2x - l) xjytZh + Xtjzh - yt)) - 2ytZh \ 
Zh&6{&6 - {x - y){xt - yt)) Q 2 {zh) /’ 


^QiUN^bh) — 


9‘^iN,TAN,bb ( yxtyt{xxj{yt - Zh) + &6Zh{x 
GAn'^mlo y ZhQ2{zh)&6 


2 )) 


+ 2 In 


xtyt 


+ 21n 


XtVtZh 1 \ 

&2{Zh)\) 


where 


0 i(a;,A) = -(-1 + y - yA)w - x(yA + w - wA) 

02(w) = (xt + y(l - xt))ytw - xxt(yyt + (yt - l)w) 
02(w) = 02(w,xt^yt) 
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03 (w) = + y) - y)ytiv + xxtiyyt + uj{-l + yt)) 

04(w) = 1 — X + XUJ 
®5("^) = X + A (1 — x) 

06 = (xt + yil - xt))yt + xxt{l - yt) 

©7 = {y{yt - 1 ) - yt)zH + x{yyt + {yt - 1 )zh) ( C . 5 ) 

08 = 2 / t ( 2x^(l + xt){yt - 1 ) + xtiy{l - yt) + yt) + x{2{l -y + yt) 

+ xt(l - 2y{l - yt) - ?,yt))) 

09 = -xj (-1 + X + y){-yt + x{2yt - l)){yt - Zh) - xtytZh{x{l + 2x) - 2y) 
+ yl{xt{x^ - y{l - x)) + (1 + x)(x - y)zh) 


and 


2 2 2 2 
mf mf mf mt 

mH± mjjo 


mr 


= 


m 


AO 


The explicit forms of the functions Fi( 2 )(yt), (^ 1 ( 2 )(yt). Hi{yt) and Li{yt) in 
Eq.(C.2) are given as 


Fi{yt) 

F2{yt) 

Gi{yt) 

G2{yt) 

Hi{yt) 

Li{yt) 


yt{7 - 5yt - 8y‘^) , y^{3yt-2) 


+ 


72{yt - 1)3 12{yt - 1)^ 


Inyt , 


yt{7>yt - 3) yt{-8yt + 2) 

12(yi-l)2+ 6(yi-l)3 ' 


yt{-yt + Syt + 2) 
24(yt - 1)3 
yt{yt - 3) 


+ 


-yt 


+ 


yt 


4(yt - 1)2 2{yt - 1)3 
1 — -isiv^dw xyt 


4(yt -1)"^ 
7a yt , 


7a yt , 


siv?6 


w 


1 


yt 


yt - 1 {yt- 1)2 
47^t^ - 79^f + 38 _ 3yl - 6t/t + 4 
108(yi - 1)3 


7a yt 


xyt 


+ 


18(yt - 1)4 

1 


Inyt 


sin^Ow 8 L yt - 1 (yt - 1)2 


Inyt 


(C.6) 


Finally, the initial values of the coefficients in the model III are 

GiHDMfjn^'^ = Gf^{mw) + Gl^{mw), 

Gf,^^{mw) = [' dx dy {C^° ((el,«)') + G§° (C^^) + G^° (/) + ((el,«)') 

^0 ^0 
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+ + CQiiCN,tt) + CQi{^N,bb))^ 

Cq^^^imw) ^ i(^N,tt)^) + C'qj (i^N,tt)'^) + C'qj iCN,tt) + C'qj {CN,bb)) 

Cl^,^^{mw) = 2^ Ow 

Cq^^^imw) = - ■ 2fl (Cq^^^(mw) - Cq^^^(mw)) 

^ m-r sm Uw 

C^q^^^(mw) = 0 , i = 5,...,10. (C.7) 

Here, we present Cq^ and Cq^ in terms of the Feynman parameters x and y sinee 
the integrated results are extremely large. Using these initial values, we ean ealeu- 
late the eoeffieients (y) and (y) at any lower seale in the effeetive 

theory with five quarks, namely u,c,d,s,b similar to the SM ease [19, 23, 25, 26]. 

The Wilson eoeffieients playing the essential role in this proeess are (fi), 

CgHDMj^^^^ (qmDM p^j. completeness, in the 

following we give their explieit expressions, 

= C^^^^{y) + Qd{Ci^^^{y)+N,Ci^^^{y)) , 

where the LO QCD eorreeted Wilson eoeffieient {y) is given by 

CLO,2HDM^^^ = ^16/23p;,2HI)M(^^)^(g/3)(^14/23_^16/23)^2HDM(^^ 

+ (C.8) 

i=l 

and y = as{mw)/o:s{y), hi and a* are the numbers whieh appear during the 
evaluation [26]. 

(y) eontains a perturbative part and a part eoming from LD effeets due to 
eonversion of the real cc into lepton pair 

Clff{ix) = Cl^^\yi) + Yresonis) , (C.9) 

where 

cT‘(m) = c|«™(M) 

+ h{z, s) (3Ci(/i) + C 2 {y) + 3C3(/i) + C 4 (/i) + 3C5(/i) + Cq{ij,)) 

- i/i(l,s)(4C3(^) + 4C4(//) + 3C5(^) + C6(M)) (C.IO) 

~ 2^^^’ + ‘^Ci{y)) + — (3(73(^) + Ci{ii) + ‘iC^{y) + Cq{ii)) , 
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and 


Y 

1 rf 


, , _ 3 ^ irTiVi ^ e+r)mY 

2 ^ / V 2 I ■ T 

"em 9 “ mVi + tmViTVi 


(3C'i(/i) + C'2(/i) + 3Cs{fi) + + 3C5(/i) + C'6(/iXP-ll) 

In eq.(C.9), the functions h{u, s) are given by 


8 mb 8 8 4 

n(u,s) = —In-inuH—-H—x 

^ ^ 9 // 9 27 9 


(C.12) 


--(2 + x)|l -X 


1/2 


In 


yi—x+i 

yr^-i 


— ITT 




2 arctan ■ 


/x-i’ 


for X = < 1 

for X = ^ > 1, 


8 8 mb 4 4 

n 0, s = —-In-In sH— ztt 

^ ^ 27 9 /r 9 9 


(C.13) 


with u = ^. The phenomenological parameter k in eq. (C.ll) is taken as 2.3. 


In Eqs. (37) and (C.II), the contributions of the coefficients C'i(/i), ...., Cq{^) are 
due to the operator mixing. 

Finally, the Wilson coefficients Cq^ (/r) and Cq^ (/i) are given by [15] 


CQiilJ) = r] CQ,{mw) , i = 1,2 . 


(C.14) 


Appendix D 

Some functions appearing in the expressions 
We parametrize the functions Ai, yl 2 , Bi, B 2 , Fi and F 2 in Eqs. (4) and (6) as 

+iAf^ A? 

A 2 = 

Bi = Ciog, 

B2 = Cio f, 

F = 2m.Cio + ^ ’ 


with 

1 ( 1 ) 




Af^ = gRe(q 


e//l 




(DA) 
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.(1) 

.(2) 

.(3) 

fW 


_p(l) 

_p(2) 


2mb 1 rf/ . 16 . 8 14 16 

=-^51-yt + -V^^)G2{yt)) , 

mb mt o 

= ^ /;yi ^ A), 

= 

= Af\gi fi), 

= 7,-12/23 + 

— 12/23 ^1 ^1—a; 

= / dx / , 

CiV,6fe A ^0 

= £dx £~''dy + + , 

— 12/23 ^1 ^1—a; 

= -T^ / / rfyC'SCelfeb) • (D.2) 
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Fig. 1. ApB as a function of sin0 for = 40mh and \rtb\ < 1 without NHB effects. 
Here ApB is restricted in the region between solid (dashed) curves for > 0 < 

0). Dashed straight line represents the SM prediction. 



Fig. 2. The same as Fig. 1, but including NHB effects with = 10 ttIt- 
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Fig. 3. Afb^^ a function of bbl'nib for = 10 mr, sin0 = 0.5 and \rth\ < 1. Here 
ApB is restricted in the region between solid (dashed) curves for > 0 < 0). 

Dashed straight line represents the SM prediction. 



nibo/mAO 


Fig. 4. The same as Fig. 3, but as a function of 
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Fig. 15. Acp{Afb) as a function of for = 40 sin6l = 0.5 and \rtb\ < 1. 
Here Acp{Apb) is restricted in the region between solid (dashed) curves for > 0 
< 0). 




Fig. 16. The same as Fig. 15, but for rtb > 1 with = 0.1 mt,. 






0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.9. 

Fig. 17. Acp{Apb) as a function of for^^^^ = IOtoi-, Cw b& = 40mf„ 

sin0 = 0.5 and \rtb\ < 1. Here Acp{Apb) is restricted in the region between solid 
(dashed) curves for > 0 < 0). 



Fig. 18. The same as Fig. 17, but for ra > 1 with = ttIt and = 0.1 mb. 








